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We report on simulations of QCD with many flavors of degenerate quarks, the DBW2 gauge ac- 
tion and naive staggered fermions, using the rational hybrid Monte Carlo algorithm. We primarily 
focus on eight degenerate quark flavors where a variety of values of the coupling constant and 
quark mass have been used in the simulations. The scaling behavior of the hadron spectrum and 
the string tension of the heavy quark potential is studied, to probe whether the zero temperature, 
C*~) continuum limit of the theory breaks chiral symmetry. 
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1. Introduction 

It is well known that vector-like gauge theories lose asymptotic freedom if the number of 
massless fermion species, Nf, exceeds some critical value Nf, where Nf = 16.5 for QCD. For Nf 
large, but smaller than Nt, a weak coupling, infrared fixed point exists in the two-loop beta func- 
tion [1, 2], which indicates that as Nf is increased from zero an interesting conformal phase may 
exist [3] before Nf becomes large enough to lose asymptotic freedom. The range of Nf where such 
a conformal phase may exist is referred to as the conformal window, and lattice simulations pro- 
vide an ideal method to determine the location of this window. Using the Schrodinger functional 
method, the running of the coupling constant can be measured on the lattice [4], which has given 
evidence for a conformal window for QCD for 12 < Nf < 16. Additional studies of the behav- 
ior of the finite temperature phase transition as a function of Nf can be used as a probe into the 
conformal window [5]. In this paper, we investigate the question of whether Nf = 8 QCD is in 
the conventional, chirally broken phase of QCD, or in the conformal phase. Previous simulations 
of eight flavor QCD [6] revealed lattice artifacts for the lattice spacings accessible at that time and 
prevented a clear statement about the phase of the theory at a single lattice spacing, and no informa- 
tion about the continuum limit. Here we use conventional, zero temperature lattice simulations to 
measure the chiral condensate and light hadron spectrum for a variety of quark masses and lattice 
spacings, so we can study the zero quark mass limit and the continuum limit. 

2. Simulations and results 

2.1 Choice of algorithm: RHMC 

We chose the Rational Hybrid Monte Carlo (RHMC) algorithm [7, 8], because it is exact 
within machine precision, its software implementation includes many optimizations, and it is ex- 
tensively used for our other 2 + 1 flavor simulations with domain wall and staggered fermions. To 
check it for our current task, we compared our results from a 4 flavor simulation, done as a 2+2 
flavor RHMC simulation (which requires the square root of the staggered fermion determinant) 
to previous results with <J> algorithm [9]. In Tab. 1, the plaquette, (yy), and masses of 4 meson 
channels are compared. 1 The plaquette values agree at the 1 a level, while the masses differ by 2 to 
3 standard deviations, making it likely that there are long autocorrelation times in the simulations 
which are not under good control. 

2.2 Choice of lattice action: Staggered fermions with DBW2 gauge 

In this exploratory work, we seek to understand the basic properties of Nf = 8 QCD at zero 
temperature, in the chiral and continuum limits, which requires simulations at many different pa- 
rameter values. We have chosen to use the naive staggered fermion action for its remnant chiral 
symmetry and its simulation speed. To help control the finite lattice spacing artifacts of staggered 
fermions, we have used the DBW2 gauge action, which produces smoother gauge fields at the lat- 
tice scale, for a given low energy physical scale, than other gauge actions. This smoothing of the 
gauge field might be expected to decrease the lattice artifacts, which we will address shortly. 

'in this paper, all numeric values of dimensionful quantities are in lattice unit unless explicitly shown. 
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Algorithm 




RHMC 


Run Length 


1000 ~ 4760 


1290 ~ 3840 


Acceptance Rate 


0.87 


0.60 


Measurement Interval 


5 


10 


Plaquette 


0.560130(14) 


0.560072(39) 


(Yy) 


0.0404(1) 


0.04105(35) 


m n 


0.3210(40) 


0.3127(30) 


«% 2 


0.3543(35) 


0.3495(41) 


m p 


0.4763(59) 


0.4610(31) 




0.4777(84) 


0.4561(37) 



Table 1: A comparison between RHMC and <t> algorithm results using naive staggered fermions and the 
Wilson gauge action for Nf = 4, m q = 0.015, /3 = 5.4. 
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Figure 1: A comparison of taste symmetry breaking. Left: compare Nf — 2 DBW2 results with quenched 
results from [10], where the chiral limit has been taken. Right: compare Nf = 2 DBW2 results with Nf = 2 
results from [11, 12] at m % jmp = 0.55. Here, except for the "DBW2" points, all labels refer to the type of 
valence quark used on a fixed dynamical ensemble. 



Additionally, using the DBW2 gauge action should help to compensate for the coarsening of 
the gauge fields that comes from adding more fermions. From asymptotic freedom, we know that 
adding more fermions to a lattice theory with a fixed low energy scale produces a larger coupling 
at the lattice cutoff, and hence coarser gauge fields. These coarse fields, in turn, make any lattice 
fermion formulation a poorer approximation to the continuum. 

2.3 Measuring taste symmetry breaking 

We have measured the taste symmetry breaking between the Goldstone pion and the local, 
non-Goldstone pion for our DBW2 simulations with naive staggered fermions and 2 dynamical fla- 
vors. In Fig. 1, our results are compared with results from improved staggered actions in quenched 
[10] and 2 flavor ("Staple+Naik" sea quarks with Symanzik improved gauge action) [11, 12] simu- 
lations. Although the comparisons are at different lattice spacings, it appears that the DBW2 gauge 
action has reduced the splittings seen with naive staggered fermions and the Wilson gauge action. 

2.4 Results with eight flavors 

We report on simulations of QCD with 8 degenerate quark flavors with three different values of 



3 



Lattice QCD with Eight Degenerate Quark Flavors 



Xiao-Yong Jin 



p 


Size 


m q 


Trajectories 


(Yy) 


m p 


ro 


aT l /G&f 


0.58 


16 3 x 32 


0.025 


1330 


~2760 


0.09973(27) 


0.812(11) 


4.39(56) 


1.73(22) 






0.015 


880^ 


- 1950 


0.06582(13) 


0.619(13) 


5.05(78) 


1.99(31) 




24 3 x 32 


0.025 


1060 


~3390 


0.100381(67) 


0.7832(30) 


4.126(96) 


1.628(38) 






0.015 


960^ 


-2930 


0.06652(11) 


0.6126(28) 


5.10(11) 


2.014(45) 


0.56 


16 3 x 32 


0.024 


970^ 


-4920 


0.13643(20) 


0.9431(38) 


3.19(18) 


1.259(70) 






0.016 


1040 


~ 3730 


0.10147(26) 


0.803(12) 


3.68(15) 


1.451(58) 




24 3 x 32 


0.024 


1010 


~3340 


0.13668(14) 


0.9693(69) 


3.120(48) 


1.231(19) 






0.016 


1040 


~3190 


0.10208(12) 


0.8085(93) 


3.793(97) 


1.497(38) 






0.008 


1000 


~2970 


0.06148(16) 


0.6022(73) 


4.716(92) 


1.861(36) 


0.54 


16 3 x 32 


0.03 
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~6220 


0.23100(20) 


1.258(17) 


2.197(52) 


0.867(21) 
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990' 


-5300 


0.19646(28) 


1.176(19) 


2.350(47) 


0.927(19) 






0.01 


1030 


~ 5520 


0.14464(37) 


0.993(14) 


2.849(51) 


1.124(20) 




24 3 x 32 


0.01 


1070 


~2860 


0.14393(39) 


1.022(17) 


2.830(48) 


1.117(19) 



Table 2: A brief summary of our simulation parameters and results. The value of a is determined from 
measuring ro and assuming ro for Nf — 8 has a physical value of 0.5fm. 

the coupling: j3 = 0.58, j8 = 0.56, /3 = 0.54; 2 different lattice sizes: 16 3 x 32 and 24 3 x 32; and 2 
or 3 different values of the quark mass for each coupling. Details can be found in Tab. 2, where the 
trajectory numbers shown are those trajectories where measurements are done. The length of each 
trajectory is 0.5 molecular dynamics time units, and measurements are done every 10 trajectories. 
The jackknife re-sampling method is used throughout our analysis and all the errors we present in 
this paper are only statistical errors. 

For all three j6 values, we start our RHMC simulations from both ordered and disordered 
configurations. In Fig. 2 we show the evolutions of (YY) f° r P = 0.58 and j8 = 0.54. We see 
that ordered and disordered starts have values of (YY) which agree after thermalization. The 
metastability for eight flavor simulations with the Wilson gauge action seen in [6] does not appear in 
our current work. In addition to the change in action, the earlier work used the inexact R algorithm, 
while here we use the exact RHMC. While we have not fully investigated what is responsible 
for the disappearance of metastability, it is very helpful to see that it is not present in the current 
simulations. 

The heavy quark potential is measured on our ensembles using the method in [13], and we fit 
the potential to the form 

V(r)=V -j + ar. (2.1) 

Fig. 3 shows the heavy quark potential for the ensemble with j3 = 0.56, m q = 0.008 on a 24 3 x 32 
lattice. The red curve is the fit to Eq. (2.1) over the range of data points denoted by triangle symbols, 
and the data points shown as diamond symbols are left out of the fit. We can clearly see confining 
behavior from the shape of the potential. Values of ro and r\ are obtained from the fit and we 
extrapolate them linearly to the zero quark mass limit as shown in Fig. 4. 

It is important to check whether Nf = 8 QCD has dynamical chiral symmetry breaking, which 
gives (Yy) a non-zero value in the zero quark mass limit. In Fig. 5, we show linear extrapolations 
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Figure 2: Evolution of (fy). Left: /3 = 0.58, m q = 0.015. Right: J3 = 0.54, m q = 0.01. In both plots, the 
data is binned in blocks of 50 trajectories. Black diamonds and red crosses represents ordered and disordered 
starts, respectively, with a lattice size of 16 3 x 32. Blue squares in the right plot are from an ordered start 
with a lattice size of24 3 x 32. 





Figure 4: Extrapolation to the chiral limit of rg and r\ . 
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Figure 6: Chiral extrapolation of pion mass. Left: Goldstone Pion mass. Right: Non-Goldstone Pion mass 
- Scalar channel which corresponds to r as ° n3 = 1 + +, of the meson propagators. 



of (yy) to m q = (left panel) and the scaling behavior of (YW){ m q = 0) versus a 2 (right panel). 
The data clearly support a non-zero chiral condensate in the continuum limit. 

If Nf = 8 QCD is in a phase with broken chiral symmetry, there should be a pseudo-Goldstone 
boson whose mass vanishes in the chiral limit. Since staggered fermions have a remnant chiral 
symmetry, we should be able to see this Goldstone particle in our simulations. We have measured 
the masses for both the Goldstone pion (jt) and the local non-Goldstone pion (712) in our simula- 
tions. In Fig. 6, the left panel shows a linear extrapolation of in\ versus m q . (Clearly there are chiral 
logarithm corrections to a simple linear extrapolation, but the coefficient of the chiral log term in 
the continuum goes as 1 /Nf, which may make the effects harder to see here than in physical QCD. 
There are also taste breaking effects in the chiral limit. We have not attempted to include any NLO 
chiral terms of this kind.) The extrapolated Goldstone mass is consistent with zero, supporting the 
conclusion that Nf = 8 QCD is in the chirally broken phase. The right panel of Fig. 6 shows the 
extrapolation of mL to the chiral limit, which is clearly non-zero for our coarsest lattice. For the 
finer lattices, the taste-symmetry breaking is much smaller. 

We have done a linear extrapolation of our measured values for m p to the chiral limit, and then 
extrapolated these to the continuum limit, where we find a non-zero value for m p (m q = 0). 
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3. Conclusions and outlook 

We have shown data from lattice simulations of QCD with 8 degenerate quark flavors that 
supports, in the chiral and continuum limits, a non-zero value for tq and r\ from the heavy quark 
potential, a non-zero value for {ijfy}, a massless % and a massive p. This argues that the conformal 
phase for QCD with Nf fermions in the fundamental representation must begin with Nf > 8. 

Having presented arguments for the behavior of 8 flavor QCD in the chiral and continuum 
limits, we are considering looking for evidence of a zero temperature conformal phase in low 
energy QCD observables with more flavors. 
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